Brownian rectifiers in the presence of temporally asymmetric unbiased forces 
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Abstract: The efficiency of energy transduction in a temporally asymmetric rocked ratchet 
is studied. Time asymmetry favours current in one direction and suppresses it in the opposite 
direction due to which large efficiency ~ 50% is readily obtained. The spatial asymmetry in the 
potential together with system inhomogeneity may help in further enhancing the efficiency. Fine tun- 
ing of system parameters considered leads to multiple current reversals even in the adiabatic regime. 
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I. INTRODUCTION 

Brownian rectifiers or ratchets are devices that convert 
nonequilibrium fluctuations into useful work in the pres- 
ence of load. Several physical models Q, 0, 0, El have 
been proposed to understand the nature of currents and 
their possible reversals with applications in nanoparticle 
separation devices 0. The possibility of enhancement 
of efficiency with which these Brownian rectifiers convert 
the nonequilibrium fluctuations into useful work has gen- 
erated much interest in this field. This, in turn, has led 
to the emergence of a separate subfield - stochastic ener- 
getics - on its own right 0, [|J . Using this formalism one 
can readily establish the compatibility between Langevin 
or Fokker-Planck formalism with the laws of thermody- 
manics thereby providing a tool to study systems far from 
equilibrium. With this framework one can calculate vari- 
ous physical quantities such as efficiency of energy trans- 
duction Q, energy dissipation (hysteresis loss), entropy 
(entropy production) jg, etc. 

Most of the studies yield low efficiencies, in the sub- 
percentage range, in various types of ratchets. This is 
due to intrinsic irreversibility associated with ratchet op- 
eration. Only fine tuning of parameters can lead to a 
large efficiency, the regime of parameters, however, being 
very narrow Recently Makhnovskii et al. [T(| con- 
structed a special type of flashing ratchet with two asym- 
metric double-well periodic-potential- states displaced by 
half a period. Such flashing ratchet models were found 
to be highly efficient with efficiency an order of magni- 
tude higher than in earlier models 0, 0, ^| . The basic 
idea behind this enhanced efficiency is that even for diffu- 
sive Brownian motion the choice of appropriate potential 
profile ensures suppression of backward motion and hence 
reduction in the accompanying dissipation. 

In the present work, we study the motion of a parti- 
cle in a rocking ratchet rocked purposefully as to favour 
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current in one direction but to suppress motion in the 
opposite direction. This is in similar spirit as in case 
of flashing ratchets proposed by Makhnovskii et al. ^(| • 
This is accomplished by applying tem porally asymmet- 
ric but unbiased periodic forcings [12, fl3l HI ] . Interest- 
ingly, such choice of forcings helps in obtaining rectified 
currents with high efficiency even for spatially symmet- 
ric periodic potentials. Still higher efficiency is obtained 
with asymmetric potentials. The range of parameters of 
operation of such ratchets is quite wide sustaining large 
loads. In addition, frictional inhomogeneity may further 
enhance the efficiency. We also see multiple current re- 
versals in the full parameter space of operation even in 
the adiabatic regime. However, multiple current reversals 
require fine tuning of the parameters. 



II. THE MODEL 

The Brownian motion of a particle in an inhomogc- 
neous medium is described by the Langevin equation [l5| 



V'(q) - F(t) k B Tj'(q) 
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where £(t) is a randomly fluctuating Gaussian thermal 
noise with zero mean and correlation, < > = 

2 5(t — t'). The periodic potential V(q) = —sin(q) — 
(/x/4) sin(2q). The parameter ^t(— 1 < fi < 1), charac- 
terises the degree of asymmetry in the potential. The 
friction coefficient 7(g) = 70 (1 — Xsin(q + 0)), with 
< A < 1 where 4> is the phase difference. F(t) is 
the externally applied periodic driving force. The cor- 
responding Fokker-Planck equation is given by 
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Since we are interested in the adiabatic limit we first 
obtain an expression for the probability current density 
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j in the presence of a constant external force Fq. The 
expression is given by 



J = 



Jo n dyl-(y) 
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where i_ (y) is given by 

i-(y) 
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It may be noted that for [i = 0, j(F ) ^ —j(—F ) for 
4> , 7r. This asymmetry ensures rectification of cur- 
rent for the rocked ratchet even in the presence of spa- 
tially symmetric potential. We assume that Fit) changes 
slow enough, i.e., its frequency is smaller than any other 
frequency related to the relaxation rate in the problem 
such that the system is in a steady state at each instant 
of time. 

We consider time asymmetric ratchets with a zero 
mean periodic driving force [T^ given by 

F(t) = \^Fo, (nr<t<riT+ir(l-e)), (5) 

= -F , (nr + -r(l - e) < t < (n+ l)r). 

Here, the parameter e signifies the temporal asymme- 
try in the periodic forcing. For this forcing in the adia- 
batic limit the time averaged current is given by 0, IT^ 
< j >= l[j(F(t)) +j(-F(t))} which is equal to 
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with 



3i = (l-e)i(l±^F ), 
h = (l + e)i(-F ). 
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The input energy £,„ per unit time is given by Q 



E in = 2 F o[(y— 
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To calculate efficiency a load L is applied against the di- 
rection of current with overall potential given by V(q) = 
—sin(q) — (/i/4) sin(2q) + qL. The current flows against 
the load as long as the load is less than the stopping force 
L s beyond which the current is in the same direction as 
that of the load. Thus in the operating range of the load 
< L < L s the Brownian particles move in the direction 
opposite to the load thereby storing energy. The average 
work done over a period is given by 
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The thermodynamic efficiency of energy transduction is 
[H, |(| r\ = E out /E in . In our present discussion all the 
physical quantities are taken in dimcnsionless units. In 
the following section we discuss the results of our calcu- 
lation. In order to evaluate currents we use the method 
of Gaussian quadrature. 



III. RESULTS AND DISCUSSIONS 

To begin with we consider a homogeneous system in 
the presence of spatially symmetric potential. In Fig. ^ 




FIG. 1: Efficiency vs load for various values of e with fixed 
F = 0.1. 




FIG. 2: Input and output energy vs load for e — 0.8 with 
fixed Fq = 0.1. The negative values of the output energy are 
not shown. The output curve is blown up ten times to scale 
with the input curve values. 



we plot efficiency as a function of load in the presence 
of temporal asymmetry e for Fq = 0.1 and T = 0.1. 
For a given e the efficiency slowly increases with load, at- 
tains a maximum and then decreases rapidly. The locus 
of peak values corresponding to different e values (with 
appropriate load L) is found to have a nonmonotonous 
behaviour with the maximum (~ 0.29) being at around 
e = 0.8 corresponding to a load of L = 0.258. For 
this value of e there is a large external driving force 
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FIG. 3: Efficiency vs load for various e with fixed Fo = 0.1. 
Inset shows the efficiency vs load for the case e = and /i = 1 . 



field, i.e., (j^)Fq(= 9Fq = 0.9) for a time duration 
t(1 — e)/2(= O.lr), in one direction which causes consid- 
erable reduction in the barrier height in that direction. 
This makes the particle cross the potential barrier easily 
and slide down the slope before the field changes. Dur- 
ing the other part of the period of F(i) (= — Fo = —0.1) 
which lasts for 0.9r the potential barrier for reverse mo- 
tion remains large thus nullifying the chances of reverse 
motion. There is an optimum value of e (~ 0.8, for the 
given parameters), however, for which considerable cur- 
rent is obtained. Beyond this value of e the duration for 
which the force remains large is so small that eventhough 
the barriers may be negligible or vanishing altogether, the 
particle may not have enough opportunity to move away 
before the force gets changed. This is the reason why the 
rectifier works most efficiently for an optimum value of e. 
For finite e current against the load is obtained for load 
L < L s . L s represents the range of load for which useful 
work is performed and is an increasing function of e. 

The useful work so obtained E out and the input energy 
Ei n are shown in Fig. [21 for a representative value of e 
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FIG. 4: Efficiency vs temperature with /i = 1 for (i) e — 
0.2, F = 0.8, (ii) e = 0.4, F = 0.5 and (iii) e = 0.8, F = 
0.1. Inset shows current as a function of temperature for 
e — 0.8, Fo = 0.1 in the absence of load. 
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FIG. 5: Efficiency vs (j> for e = 0.8, fi = 1,L = 0.44 for (i) 
A = 0.0 and (ii) A = 0.9 with fixed F = 0.1 and T = 0.1 
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FIG. 6: Current vs temperature curve showing two current 
reversals for /i = —1, A = 0.9 and e = 0.34 with <f> = 
1.005-7T, Fo — 0.3 and L — 0. Inset shows the current in the 
presence of lone asymmetry parameters (A, e, p). 



as a function of load. The input energy decreases to 
a minimum value for a load larger than L s . Moreover, 
it remains positive, as expected, over the entire range. 
The output energy shows a peak with load in the region 
where the input energy is monotonously decreasing. It 
then becomes negative for L > L s as anticipated. The 
qualitative behaviours of efficiency and energies shown in 
Figs. and [5] are similar to those in reference 01 for 
the flashing ratchet. 

In Fig. 01 we consider the periodic potential V(q) to be 
spatially asymmetric together with a temporally asym- 
metric external driving force field. The potential asym- 
metry enhances the efficiency of energy transduction as 
well as widens the range of load. This is due to the fact 
that for f / the presence of asymmetric parameter 
//(> 0) further reduces the potential barrier for forward 
motion and enhances the barrier for backward motion. 
Moreover, as can be seen from the inset, one can get fi- 
nite current even when e = with finite stopping force L s 
in contrast to the symmetric potential case. From Figs. 
Hand |2I it is clear that the temporally asymmetric forces 
not only enhance the efficiency of energy transduction 
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but also widen the operation range of load against which 
the ratchet system works. 

In Fig.^Jwe plot efficiency as a function of T for various 
e values in the presence of potential asymmetry (ji > 0). 
The efficiency decreases with temperature. The rele- 
vant physical parameters chosen for optimal efficiency 
are mentioned in the caption. From the inset it is to be 
noted that the current peaks as a function of temperature 
yet efficiency decreases monotonically. This implies that 
thermal fluctuation do not favour energy transduction in 
this case. It is worth mentioning that for given temper- 
ature and e the efficiency shows peaking behaviour as a 
function of Fq ; the efficiency being zero for Fq — as well 
as for large Fq for in these limits output current vanishes 
in the absence of load. 

Next, we present the effect of frictional inhomogeneity 
(7 = 7(9); A 0). In Fig. we plot the efficiency as a 
function of the phase difference between the potential and 
the friction coefficient 7(5) for a typical case. We observe 
that the inclusion of this parameter A further increases 
the efficiency in a range of 4> depending on other parame- 
ter values. It is worth mentioning that for inhomogeneous 
systems the efficiency peaks with temperature in a lim- 
ited range of parameters. With frictional inhomogeneity 
the range of temperature in which one can obtain output 
current with finite efficiency is extended to a large tem- 
perature where contribution of A dominates over other 
parameters. In Fig [fj] we show that by properly choosing 
the parameters we can obtain multiple current reversals 
as a function of temperature. It should be noted that 
such reversals are not possible in the homogeneous case 



in the adiabatic regime |17| . The inset shows current as 
a function of individual parameters (e , fi , A). The plots 
indicate that individual parameters cannot bring about 
current reversals separately. However, the possibility of 
current reversals arises due to the combined effect of the 
three asymmetry parameters considered. We have also 
observed more number of current reversals than shown 
in Fig. by fine tuning the parameters. 



IV. CONCLUSIONS 

We find large efficiency for rocking ratchets driven 
by temporally asymmetric periodic field the origin of 
which can be traced to the suppression of backward 
motion. The observed efficiency is much higher than the 
earlier reported values eventhough the ratchet operates 
in an intrinsically irreversible domain. This asymmetry 
factor has also helped in increasing the range of load 
of operation of the ratchet. We also observe multiple 
current reversals in the adiabatic limit by proper fine 
tuning of different parameters. These reversals are 
attributed to inherent complex dynamics of the system. 
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